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Phase-Field Reaction-Pathway Kinetics of Martensitic Transformations in a Model Fe 3 Ni Alloy 
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A three-dimensional phase-field approach to martensitic transformations that uses reaction pathways in place 
of a Landau potential is introduced and applied to a model of FeaNi. Pathway branching involves an un- 
bounded set of variants through duplication and rotations by the rotation point groups of the austenite and 
martensite phases. Path properties, including potential energy and elastic tensors, are calibrated by molecular 
statics. Acoustic waves are dealt with via a splitting technique between elastic and dissipative behaviors in a 
large-deformation framework. The sole free parameter of the model is the damping coefficient associated to 
transformations, tuned by comparisons with molecular dynamics simulations. Good quantitative agreement is 
then obtained between both methods. 

PACS numbers: 64.60.-i, 64.70.K-, 81.30.Kf 



(N 



CO 



Cd 



T3 

C 
O 

o 



> 

in 

in 

o 
o 



X 



Nanoscale materials that undergo martensitic transforma- 
tions (MT) bear the promise of an exceptional technological 
revolution [1]. MTs are displacive structural transitions as- 
sociated to large inelastic strains, that occur under tempera- 
ture or loading changes, from a high- (austenite) to a low- 
symmetry state (martensite) declined in a number of "vari- 
ants" |2I, of time scale down to subnanosecond order [|3|,|4[]. 
Bulk MTs in large samples lead to complex microstructures, 
due to competing long-range elasticity, and crystallographic 
constraints on variants yfl. Kinetics of MTs has been investi- 
gated at small scales by molecular dynamics (MD) (e.g., J5|]), 
whereas continuum-mechanics-based phase-field (PF) models 
IsLZI] must be used for large sizes and simulation durations. 

The unsolved issue addressed in this Letter consists in seek- 
ing quantitative agreement between PF, and MD in its opera- 
tive range of size- and time-scales, in a time-dependent set- 
ting. We focus on the illustrative case of strain-driven trans- 
formations near K in a stoichiometric (ordered) FeaNi alloy 
stable at low temperatures only 18l|9[], that undergoes a proper 
(i.e., with no shuffling) austenite 7(fcc) — > martensite a(bcc) 
transformation along a path of homogeneous deformation of 
the unit cell. Consistently benchmarking PF calculations by 
MD simulations requires adjusting the PF model using the em- 
pirical potential of the simulations, instead of more accurate 
first-principles methods (e.g., BlfJIl ). Disregarding magnetic 
degrees of freedom, we use a Meyer-Entel EAM potential de- 
veloped to investigate the phase diagram of the MT transition 
in Fe^Nii-a; alloys yllsfl (but see also GlJ])- 

In the PF method for proper MTs, the nonrelaxed Helm- 
holtz energy density has been modeled by a Landau poten- 
tial for the total strain. Levitas et al. recently extended this 
formulation to large strains, using a vector order parameter r\ 
associated to a Landau potential that describes the transforma- 
tional part of the strain |7|]. However, due to group-subgroup 
relations in the lattice symmetry point group (PG), the fee 
— > bec transformation is reconstructive 11111 . That is, once 
a martensite variant is reached from the parent phase, differ- 
ent austenite variants, among which the original one, can be 
reached in turn from the martensite (Fig. H}. Repeatedly ap- 
plying PG transformations thus implies considering an infinite 



set of variants. Landau theory is then inapplicable, though an 
approximate theory with non-commensurable order parameter 
can be used 11111 . 

In this context, as a simple alternative to using Landau po- 
tentials, we introduce a PF approach based on reaction path- 
ways (PF-RP). The RP is a minimum-energy path that links 
two (meta)stable states via a saddle point (e.g., 11211 ). Con- 
sider two states, austenite (A) and martensite (M), of defor- 
mation gradients F = I + Vu = F A ' M (I is the identity and u 
is the material displacement) with respect to an arbitrary refer- 
ence state, which are local energy minimizers with associated 
elastic moduli tensors C A ' M , computed by molecular statics 
(MS). Although the actual path might slightly differ (121, we 
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FIG. 1. (color online). Left: cycling transformations austenite A 
— >• martensite M — s> A, etc., produces in strain gradient (F) space 
an infinite backbone of preferred RPs for the transformational strain 
Ft. However, Ft can depart from it during out-of-path transitions. 
Right: Inelastic energy along the line (dash) that minimizes the dis- 
tance between two nonconnected pathways, in a hypothetic pathway 
arrangement chosen for ease of representation (top). The role of a 
and 77 are emphasized, inelastic energies on pathways / in (%), con- 
stant in this case, being set to 0. The parameter a scales the barrier 
energy, taken proportional to the distance between RPs, and n is an 
empirical shape parameter (bottom), see Eqs. 10 and 10). Thus, di- 
rect transitions between remote RPs are inhibited. 

approximate the RP between these states by optimizing over 
volume 0] along a Bain path 0] F B (s) = sF A + (1 - s)F M , 
where < s < 1 is a path coordinate. Thus along the RP, the 
Helmholtz energy density is f(s) = min^n / MS (fcF B (s)) 



where f MS is the energy density from MS, associated to the 
deformation gradient F(s) = fc(,s)F B (s) where k(s) is the 
volume-minimizing factor. To comply with crystal symme- 
try, this first RP is duplicated and rotated for any rotations R 
of the austenite and martensite point groups (PG). Because 
of the high symmetry of the considered phases, only 3 of the 
24 possible rotations lead to new RPs (e.g., rotations R, 



ted by MS) 
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of the first martensites, |2|[|). As a substitute for the Landau po- 
tential of the total strain, we gather these RPs in a backbone 
construct for the potential. Limiting ourselves to moderate 
deformations (up to 50%), we consider only the first 21 vari- 
ants 111 311 . The reference state is the austenite denoted by fee 
where F A = I. This backbone captures the most important 
information about energy barriers. 

Quite generally, f(s) embodies the elastic energy near its 
minima s = 0, 1, so that s is not exclusively related to trans- 
formational behavior. To deal with a large-deformation the- 
ory where Landau-like parameters are associated to a trans- 
formational strain |7||, we adapt a small-strain procedure (in- 
troduced in 01411 and discussed in 111511 ) allowing one to derive 
a transformation-related potential that excludes linear-elastic 
energy. An internal transformation variable rj (the phase- 
field) is introduced by^ defining the transformational gradient 
along the RP to be F t (i]) = F(s(rj)), where the function 
r](s) that provides by inversion the relationship s = s(r]) in 
this definition must be determined. Adopting the usual large- 
deformation multiplicative composition law, define the elastic 
gradient as F e (s, rf) = F(s) • F^ (77), and write the energy 
density along the RP in the alternative form 

f(s,r)) = iE e ( S ,77) : C(«) : E e (s, V ) + f m (ri). (1) 

The first term is an elastic energy expressed using the Green- 
Lagrange elastic strain E e = | (Fj • F e — I) , and /„, is the 
inelastic energy of the RP 114n . To avoid brutal elastic vari- 
ations, the elastic tensor C is made s-dependant along the 
path and taken as a cubic interpolation between C A ' M with 
C'(s) = at both ends [7!]. Functions fi n (r]) and 77(5) are 
then obtained by imposing an exact equality between /(s) 
and /(sj) in a relaxed state where rj is adiabatically elim- 
inated (153, which yields the necessary equations to be solved 
numerically: 



f( S ) = f(s, V ), d v f(s,r,) = 0. 



(2) 



Similarly, a mechanism whereby a phase strain can leave a RP 
to rejoin a neighboring one has to be added [J. A transfor- 
mation gradient F t outside RPs is introduced as a generalized 
internal variable and is associated to a new potential. Keeping 
in mind that this potential corresponds to transient states be- 
tween two arbitrary RPs, (i.e., f in should be a function of F t ) 
the energy to be added to f> (%) is simply chosen as pro- 
portional to the distance rffc(Ft) = min^ |F t • F t {rjk) — I| 
from the k lh RP, where |A| = {A^A^fl 2 (Fig.[T] left). The 
contribution for one RP introduces one parameter a (to be fit- 
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where r/j, is the argmin in dk (F t ), a rotation-independent (i.e., 
objective) function. 

For the complete pathway tree, the overall inelastic energy 
is an interpolation between potentials 



/ in (F) = ^ Wfc (F t )/ i ( 1 f ) (F t ) 



(4) 



with a partition of unity u>k chosen so that a RP dominates 
its immediate surrounding, i.e., Wk = 1 for F t near the (fc)th 
RP. A simple and convenient choice is to use a function of the 
distance dk defined above: u>k = d" k n /J2i d^ n with n > 
controling the transition between pathways. Whereas n — > 00 
makes f m switch to the nearest RP, best agreement with MS 
is obtained using n m 2, which provides smoother transitions 
(Fig.[U right; see also 11611 ). The full potential now reads 

/(F,F t ) - |E e (F,F t ) : C(F) : E e (F, F t ) + / in (F t ), (5) 

where C(F) is interpolated from the C(sfe)s, using an equa- 
tion similar to Eq. (0), with d k {F) — min s& F-F^ ~ 1 (s fc )- 
I| to define s^. Along RP k, we have Wk = 1 and w%^k = 0, 
and Eq. (O is equivalent to Eq. ([TJ. For transition between 
pathways, the energy barrier is proportional to the distance 
between RPs, which naturally inhibits unphysical transitions 
between "distant" variants (Fig. Q] right). This approach is 
quite different from the interpolation scheme used in 121] and, 
we believe, simpler to handle, at least for reconstructive trans- 
formations involving an extended reaction tree. 

By construction, linear-elastic energy is removed from 
/i n (F t ) which, apart from small nonlinear elastic contribu- 
tions, describes the non-convex (unstable) part of the energy 
along the RP 11511 . Hence, damping can be prescribed for F t 
while leaving linear-elastic wave dynamics undamped, consis- 
tently with the negligible character of viscoelasticity in solid 
metals. With / given by (O, F t follows by hypothesis a time- 
dependent Ginzburg-Landau kinetics of parameter v. 



-^/(F.Ft). 



(6) 



In v are lumped dissipation mechanisms such as vibrational or 
magnetism entropy 11711 . nonlinear acoustic waves, as well as 
couplings to inessential lattice degrees of freedom that were 
adiabatically eliminated when computing f(s). 

Finally, the dynamics of u obeys the equation pii = V • cr, 
where p is the local density, and where the Cauchy stress a is 
related to the first Piola-Kirchhoff stress P = <9/(F, F t )/<9F. 
The model is implemented in a Lagrangian code using an 
element-free Galerkin (EFG) formulation 11811 in total strain, 
with explicit time integration in a form able to handle acoustic 
wave propagation and rapid phase changes 11911 . This least- 
square formulation of EFG produces smooth fields for u and 
F t with no pinning at interpolation nodes. Hence, including an 



interface-penalizing (gradient) term in the energy is not neces- 
sary, the same overall effect being obtained by keeping finite 
the distance between interpolation nodes. Contrary to a and 
n, the free parameter v must be fitted on global simulations. 

A benchmark MD simulation is conducted with initial tem- 
perature T — K. An austenitic cube of size L = 30ao 
with lattice parameter ao = 3.64 A (108 000 atoms) is de- 
formed with time t < if, with if = 80 ps, according to a 
time-dependent overall strain F(i) = I + (t/tf) [aF Ml + (1 - 
a )F M2 -I] with components of F Ml = (1.105; 1.105; 0.789) 
and F M2 = (0.789; 1.105; 1.105) and a ranging from 0.1 to 
0.3 to obtain various compounds of two martensite variants 
Mi 2 in the final configuration. Periodic boundary conditions 
(PBCs) are used. At i = if, acoustic waves have run 26 times 
across the sample, which evenly spreads out perturbations. 

For the corresponding PF-RP simulations, node spacing 
controls the spatial resolution. A good compromise between 
resolution and computational cost is obtained with 27 atoms 
per interpolation node (4000 nodes) placed in the same fee 
arrangement as that used for MD. Values a = 1.9 GPa and 
n = 2.2 lead to a good reproduction of the energy between 
two martensite variants M\ and M.%, for prescribed strains 
F = qF Mi + (1 - a)F M2 with a ranging between and 
0.5 (see d). 

For PF-RP and MD, F is used to monitor deformation, 
and to identify variants. For MD it is obtained from dis- 
placements u of neighboring atoms, by least-square optimiza- 
tion. Final states in MD are made of bands of martensite for 
a < 0.2 whereas "chessboard"-like structures (20] emerge for 
a = 0.3. For PF-RP, the viscosity v controls the final states, 
a very good match between MD and PF-RP being obtained 
with v < 14 mPa.s (Fig.|2]i. Lowering v does not change the 
final states, but increasing it turns the bandlike structures into 
chessboard (0.014 < v < 0.06), then into a homogeneously 
deformed state (v > 0.06). Phase volume fraction are moni- 
tored, see FigOJb)- The ti me to nucleation (TTN) is delayed 
by increasing v, best match between MD and PF-RP TTNs 
being obtained with again v = 14 mPa.s. 

This value of v is used in all calculations below. It is 
within a factor ~ 2.5 of Fe and Ni viscosities at their melt- 
ing point, close to one another (5.8 mPa.s and 5.4 mPa.s, re- 
spectively 112111 ). An explanation may be that configurational 
changes associated to barrier-crossing involve large atom mo- 
tions comparable to that encountered on the liquidus (even 
though martensitic transformations are non-diffusive), with 
similar damping effects. 

Analogous chessboard patterns have been observed and re- 
produced by PF in connection with diffusive alloy decompo- 
sition 02211 . For displacive transformations, these ubiquitous 
structures Il20ll have also been obtained in two-dimensional 
PF calculations: very unstable, they are stabilized in small 
samples, but decay into more conventional laminatelike struc- 
tures at large sizes 12311 . We make here a first exploration of 
this physically important effect in three dimensions (3D) us- 
ing PF-RP for an imposed deformation F(i) with a = 0.3. 
Because of inertial dynamics, convergence of strains to stable 
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FIG. 2. (color online). Comparison of deformation measure between 
MD (top) and PF-RP (bottom) for a = {0.1, 0.2, 0.3} at 80 ps (light: 
martensite; dark: incompletely-transformed austenite). 



values is limited by wave propagation. Therefore, the ratio 
tf/L between simulation time if and sample size L must be 
kept constant to allow for meaningful comparisons. Taking 
the size and duration of the previous calculation as a refer- 
ence (j3 = 1), L and if are increased by factors f3 = 2 to 
/? = 7. As deformation proceeds and for /3 > 3, the initial 
austenite goes through a chessboard state that decays, via an 
intermediate mixed structure, to a complex three-dimensional 
laminate state of austenite mixed with twin bands of marten- 
site. This sequence is illustrated on Figs. |2a) for the largest 
size L = 210 ao (/3 = 7) and time if = 560ps; see IU6I1 
for animated sequences. Interestingly, for (3 = 7 variant 3 
is produced at intermediate times (Fig. [5J)) and vanishes at 
the end of the simulation. Indeed, the prescribed strain F(i) 
is defined as an average of initial austenite [strain I, volume 
fraction (VF) 1 — t/tf], and of the two variants [strain F 1 , 
VF at/tf and F M2 , VF (1 - a)t/t f ]. However, the impor- 
tant energy gain when martensite forms (-13 meV/atom) fa- 
vors larger fractions of both variants 1 and 2. Noting that the 
average strain produced by a combination of the three variants 
is null, this is balanced by a "back" strain proportional to F 3 
inducing the formation of the third variant. 



In the final state, the martensite compound forms a two- 
dimensional structure more complex than a simple laminate, 
in which the two possible orientations of twin interfaces con- 
sistent with boundary conditions at habitat planes 12[] are si- 
multaneously present. At meeting points of 90° -related inter- 
faces, these boundary conditions cannot be satisfied and high 
elastic strains result. Relaxation occurs through a moderate 
formation of "reversion" austenite (less than 0.05%) of the 
Mi — > A2 path, see Fig.[TJ represented as (barely noticeable) 
dark regions in Fig. [3ja), right. Additional MD simulations 
and PF-RP calculations were made, using PBCs with planes 
rotated by an angle of 5°. This small change inhibits chess- 
board patterning whatever the system size. A laminatelike 
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FIG. 3. (color online), (a) Intermediate and final microstructures for 
/3 — 7. (b) Volume fraction of variants for small (8 = 1) and large 
(J3 = 7) simulations for MD (dots) and PF-RP (plain and dashed 
lines). Simple composition rule between initial austenite and vari- 
ants is plotted (black lines). For /3 = 7, variant 3 appears as an 
intermediate phase. A small fraction of "reversion" austenite (label 
A2) is produced at variant angle points. 



structure takes place almost instantaneously without giving 
rise to any remarkable intermediate state. This suggests that 
stable chessboards may be difficult to observe experimentally 
for displacive transformations, even in small samples. 

To conclude, we introduced a dynamic phase-field tech- 
nique for martensitic transformations, fully compliant with 
crystal symmetries, that alleviates the need for vector Landau 
parameters, and obviously adaptative to a variety of situations. 
We illustrated it by an application to a model alloy. The good 
results obtained, which contrast with the crude approxima- 
tions involved in modeling the energy landscape outside RPs, 
show that the details of these "outer" regions are most likely 
inessential to the main picture and confirm the relevance of 
a reaction-pathway approach to these questions. This view is 
supported by the agreement found with MD in size and time 
domains where both techniques could be compared, with a 
gain of two orders of magnitude in computational cost in fa- 
vor of PF-RP. 

CD. and Y-P.P. thank G. Zerah and L. Truskinovski for 
seminal discussions. 
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